S=eye(5)
whos
nnz(S) 
nonzeros(S)
nzmax(S)
spy(S)
[i,j,s]= find(S)
[i,j]= find(S)
